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INTRODUCTION 


Recent  proposals  for  orbiting  space  systems  have  stimulated  research 
in  the  design,  analysis  and  control  of  large  space  structures.  Such 
structures  are  envisioned  to  be  very  large  (one  review  paper  [1]  speculates 
about  "a  structure  nearly  the  size  of  Manhattan  Island")  and  are  often  in 
the  form  of  lattice  structures.  A  lattice  structure  is  a  network  of 
slender  elements  which  are  connected  by  joints.  This  paper  is  concerned 
with  the  analysis  of  the  dynamics  and,  in  particular,  the  wave-like 
motion  of  lattice  structures. 

It  is  convenient  to  use  matrix  methods  to  analyze  the  dynamics  of 
lattice  structures.  In  particular,  the  dynamics  of  a  lattice  member  may 
be  described  by  a  transfer  matrix  relationship,  and  the  dynamics  of  a 
joint  in  a  lattice  structure  may  be  described  by  a  joint  coupling  matrix 
relationship.  In  order  to  study  wave  propagation  in  lattice  structures, 
it  is  convenient  to  introduce  wave-mode  coordinates  and  scattering  matrices 
Wave-mode  coordinates  are  obtained  through  the  process  of  diagonalizing 
or  uncoupling  the  transfer  matrix  relationship.  When  the  transformation 
used  to  uncouple  the  transfer  matrix  is  applied  to  the  joint  coupling 
matrix  relationship,  the  scattering  matrix  is  obtained. 

In  this  paper,  the  derivations  leading  to  wave-mode  coordinates  and 
scattering  matrices  in  lattice  structures  are  reviewed.  In  addition, 
simple  one-dimensional  examples  are  given  to  illustrate  how  wave-mode 
coordinates  and  scattering  matrices  may  be  used  to  describe  dynamics  and 
wave  propagation  in  large  space  lattice  structures. 


ANALYSIS 


State  Vectors  and  Transfer  Matrices 


In  the  matrix  formulation  of  the  dynamics  of  lattice  structures,  the 
state  or  configuration  of  the  lattice  is  determined  by  state  vectors 
which  exist  at  each  point  in  the  lattice.  The  state  vector  at  a  particular 
point  is  a  column  matrix  whose  components  may  be  displacements,  rotations, 
internal  forces  and  internal  moments.  The  actual  components  contained 
in  a  state  vector  at  a  particular  point  depend  on  the  model  used  to 
describe  the  latt  ce  at  that  point. 

If  the  members  of  a  lattice  structure  are  slender,  they  may  be  modeled 
as  one-dimensional  continua.  The  components  of  the  state  vectors  (and 
the  state  vectors  themselves)  of  a  one-dimensional  continuous  member  are 
functions  of  time  and  one  local  spatial  coordinate,  say,  x  which  extends 
along  the  length  of  the  member.  At  a  particular  location  x  =  x^,  the 
Fourier  transform  of  the  state  vector  of  a  one-dimensional  member  is  given 
by  [2] 


CO 

z^x^w)  =  J  £(x1,t)e  1L 


where  £(x^,t)  is  the  state  vector  at  x  =  x^,  z_(x^,w)  is  the  Fourier  trans¬ 
form  of  the  state  vector  at  x  =  x^,  t  is  time,  to  is  radian  frequency,  and 
i  =  J-i  .  The  state  vector  £(x^,t)  is  given  in  terms  of  the  transformed 
state  vector  £(x^,uj)  by  the  inverse  Fourier  transform  relationship  [2] 


00 

z.Cxj.t)  =  J  _z(x1,w)e1L 
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Throughout  this  paper,  an  overbar  will  signify  a  Fourier  transform,  and 


an  underbar  will  signify  a  vector  or  a  matrix. 

The  Fourier  transforms  of  the  state  vectors  at  two  locations  x  =  x^ 
and  x  =  x^  of  a  one -dimensional  lattice  member  can  be  related  by  a  transfer 
matrix  according  to 

z(x2,u))  =  T(x2  -  x1,w)z_(x1,u>)  (3) 

where  T(x2  -  x.  ,u>)  is  the  n  x  n  transfer  matrix  (here  n  is  the  (even) 
number  of  components  of  the  transformed  state  vectors  z^x2,w)  and  z^x^.co)), 
which  is  a  function  of  the  separation  (x2  -  x^)  and  the  frequency  to.  The 
transfer  matrix  is  derived  from  the  equations  of  motion  and  the  constitutive 
equations  of  the  member.  Transfer  matrices  for  several  specific  types  of 
lattice  members  are  given  in  [3] . 

Wave-Mode  Coordinates  and  Propagation  Constants 

The  transfer  matrix  T_  in  eqn.  (3)  generally  contains  nonzero  elements 
off  its  main  diagonal,  and  therefore  the  scalar  equations  relating  the 
components  of  z_(x2>uj)  and  z^(x^,oa)  are  coupled.  In  order  to  study  wave 
propagation  in  lattice  structures,  it  is  convenient  to  cast  the  transfer 
matrix  relationship  given  by  eqn.  (3)  into  a  different  form.  In  particular, 
it  is  desirable  to  diagonalize  the  transfer  matrix,  and  thus  produce  an 
uncoupled  set  of  scalar  equations.  This  diagonalization  may  be  accomplished 
by  the  following  procedure. 

First,  the  eigenvectors  of  T  are  found,  and  the  n  x  n  wave-mode  matrix 
Y(o.)  is  assembled.  The  columns  of  Y(ui)  are  the  eigenvectors  of  the 
transfer  matrix  T.  Next,  a  new  vector  w(x,u)  is  defined  by  the  equation 
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w(x,to)  =  Y  (to)  z^x.to) 


The  vector  w(x,w)  will  be  discussed  and  interpreted  shortly.  From  eqn.  (4) 
it  follows  that 


z(x,to)  =  Y(to)w(x,w) 


Note  that  the  existence  of  Y  (oo)  in  eqn.  (4)  depends  on  the  linear 
independence  of  the  eigenvectors  of  T  [4],  Substitution  of  eqn.  (5)  into 
eqn.  (3)  gives 


Y(u)w(x2,cl>)  =  T/x2  -  x1  ,w)Y(cu)w(x1  ,w) 


Premultiplication  of  both  sides  of  eqn.  (6)  by  Y  (to)  gives 


w(x2,io)  =  W((jj)w(x1,w) 


where 


W(co)  =  Y  (oj)T(x2  -  x1,oo)Y(aO 


The  n  x  n  matrix  W(to)  is  called  the  wave-mode  propagation  matrix.  Assuming 
that  the  eigenvectors  of  T  are  linearly  independent  (and  that,  therefore, 

Y  1(to)  exists),  the  matrix  W(oo)  is  a  diagonal  matrix  of  the  form  [4] 


I 


W(u>) 


K 


1 


0 


0 


.  0 

0  ...OK 


(9) 


where  the  (i  =  1,2,  ...,  n)  are  the  n  eigenvalues  of  T.  For  uniform 
lattice  members,  which  are  considered  exclusively  in  this  paper,  the 
transfer  matrix  is  cross-symmetric,  and  therefore  the  eigenvalues 
of  T  occur  in  pairs  (K. ,  1 /K^ ) ,  j  =  1,2,  ....  n/2  [5]. 

The  vector  w(x,  w)  given  by  eqn.  (4)  is  called  the  wave-mode  vector, 
and  the  components  of  w(x,  oj)  are  called  wave-mode  coordinates.  The 
wave-mode  coordinates  are  given  in  terms  of  the  physical  components  of 
the  state  vector  by  eqn.  (4) ,  and  the  physical  components  of  the  state 
vector  are  given  in  terms  of  the  wave-mode  coordinates  by  eqn.  (5).  In 
general,  each  wave-mode  coordinate  will  be  a  combination  of  all  the 
components  of  the  state  vector.  The  wave-mode  matrix  Y(co) ,  which  gives 
the  transformation  between  the  wave-mode  vector  and  the  state  vector, 
is  similar  in  some  respects  to  the  modal  matrix  used  in  the  analysis 
of  multiple  degree  of  freedom  vibratory  systems. 

Because  the  matrix  W(u))  in  eqn.  (7)  is  diagonal,  the  wave-mode 
coordinates  are  uncoupled  from  each  other  as  they  propagate  from  x  =  x^ 
to  x  =  along  the  lattice  member.  The  components  of  w(x,  u')  are 
called  wave-mode  coordinates  because  they  represent,  for  many  lattice 
member  models,  propagating  waves.  Eqn.  (9)  can  be  written  in  the  form 


7 


W  (co)  = 


Y1(x2-x1) 


Y2(x2-Xi) 


.  0  e 


Yn<*2-Xl> 


O 


where  the  quantities  Yi  (i  =  1,2,  ....  n)  satisfy  the  equations 


K.  =  e 
1 


Yi(x2-x1) 


(1 


or 


ln(K  ) 


'i  (x2-Xl)  X“V,V 


(1 


The  quantities  Y^  are  called  "propagation  constants",  and  are  in  general 
complex.  Since  the  eigenvalues  occur  in  pairs  (K^ ,  1/K.),  the  propa¬ 
gation  constants  occur  in  pairs  (y^ ,  ~Yj ) ,  j  =1,2,  n/2.  The 

general  form  of  the  propagation  constant  pair  (y . ,  -y.)  is  [6] 


Yj  =  a  (to)  +  ik(u)) 
-Yj  =  -a(io)  -  ik(co) 


(1 


(1‘ 


where  ci(w)  is  called  the  attenuation  and  k(w)  is  called  the  wave  number. 

Propagating  waves  represented  by  wave-mode  coordinates  may  be 
classified  according  to  the  nature  of  their  corresponding  propagation 


constants.  If  a  propagation  constant  is  purely  imaginary,  the  wave 
represented  by  the  corresponding  wave-mode  coordinate  is  not  attenuated 


as  it  propagates.  If  a  propagation  constant  has  a  nonzero  real  part  and 
a  nonzero  imaginary  part,  the  wave  represented  bv  the  corresponding  wave¬ 
mode  coordinate  is  attenuated  as  it  propagates.  A  propagation  constant 
which  is  purely  real  corresponds  to  a  wave-mod--'  coordinate  which  does  not 
represent  a  wave  at  all,  but  rather  represents  a  nonpropagating  spatially 
attenuated  vibration.  If  the  imaginary  part  of  a  propagation  constant 
is  linearly  proportional  to  frequency  co,  the  wave  represented  by  the 
corresponding  wave-mode  coordinate  is  nondispersive,  and  the  phase 
velocity  of  the  wave  is  given  by  the  inverse  of  the  constant  of  pro¬ 
portionality.  If  the  imaginary  part  of  a  propagation  constant  is  not 
linearly  proportional  to  frequency,  the  wave  represented  by  the  corre¬ 
sponding  wave-mode  coordinate  is  dispersive. 

The  physical  interpretation  of  the  fact  that  propagation  constants 
for  uniform  members  occur  in  pairs  (y  ,  -y^ )  is  that  for  uniform  members 
identical  waves  may  propagate  in  either  the  direction  of  increasing  x 
or  the  direction  of  decreasing  x.  Propagation  constants  with  a  nonpositive 
real  part  and  a  negative  imaginary  part  correspond  to  wave-mode  coordi¬ 
nates  which  represent  waves  which  propagate  in  the  direction  of  increasing 
x,  and  propagation  constants  with  a  nonnegative  real  part  and  a  positive 
imaginary  part  correspond  to  wave-mode  coordinates  which  represent  waves 
which  propagate  in  the  direction  of  decreasing  x.  For  passive  lattice 
members  (that  is,  for  lattice  members  with  no  external  energy  inputs), 
propagation  constants  with  a  positive  imaginary  part  must  also  have  a 
nonnegative  real  part,  and  propagation  constants  with  a  negative 
imaginary  part  must  also  have  a  nonpositive  real  part  [6] . 


Joint  Coupling  Matrices  and  Scattering  Matrices 


The  dynamics  of  a  joint  in  a  lattice  structure  may  be  described 
by  a  joint  coupling  relationship  of  the  form  [7] 


=  f 
—ext 


where  the  z.  (i  =  1,2,  m)  are  the  state  vectors  of  the  m  members 

—i 

which  are  connected  to  the  joint  at  the  respective  points  of  contact 
between  the  members  and  the  joint,  ^ext  is  a  vector  containing  the 
Fourier  transforms  of  the  external  forces  and  moments  which  are  applied 
to  the  joint,  and  B  is  a  matrix  called  the  joint  coupling  matrix.  The 
joint  coupling  matrix  is  derived  from  the  equations  of  motion  of  the 
joint,  the  constitutive  equations  of  the  joint,  and  the  geometric  com¬ 
patibility  requirements  at  the  joint.  The  explicit  form  of  the  joint 
coupling  matrix  for  arbitrary  two  and  three-dimensional  rigid  joints 
with  mass  is  given  in  [7]. 

Eqn.  (15)  is  written  in  terms  of  the  state  vectors  z_  (i  =  1,2,  ..., 
In  order  to  describe  the  dynamics  of  a  joint  in  terms  of  wave-mode  coordi 
nates,  the  following  transformations  are  introduced: 


z  ,  =  Y  .  (u)w.  ,  i  =  1,2 . m 


th 


where  Y  .  (w)  is  the  wave-mode  matrix  of  the  i  member  which  is  attached 


— 1 


.  th 


to  the  ioint,  and  w  is  the  wave-mode  vector  of  the  i  member  at  the 

j  — i 

.  th 


point  of  contact  between  the  i  member  and  the  joint.  Substitution 
of  eqn.  (16)  into  eqn.  (15)  gives 


(  Y^  (o')  w^ 


Y9(u)w2 


B  < 


Y  (u)w 
— m  — m  J 


>=  f 


-ext 


which  can  be  written  in  the  form 


r  w_^  ^ 


~2 


w 

^  -mj 


■^ext 


where 


C  -  B 


'Y^co)  0 


0  Y2(ic) 


(17) 


(18) 


0  Y  (w) 
— m 


(19) 
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Eqn.  (18)  describes  the  dynamics  of  a  joint  in  terms  of  wave-mode  vectors 


w..  As  discussed  above,  the  wave-mode  vectors  contain  wave-mode  coordi- 
— 1 

nates  which  represent  waves  travelling  along  the  members.  In  describing 
the  dynamics  of  a  joint  in  terms  of  wave-mode  coordinates,  it  is  con¬ 
venient  to  group  the  wave-mode  coordinates  into  two  groups:  wave-mode 
coordinate  which  represent  waves  "entering"  the  joint,  and  wave-mode 
coordinates  which  represent  waves  "leaving"  the  joint.  A  wave-mode 
coordinate  in  the  i*"  1  member  attached  to  the  joint  represents  a  wave 
entering  the  joint  if 


•gn  I  Im(‘) ) 


(n  . 

—  i 


e 

“XI 


(20 


where  ')  is  the  propagation  constant  cor  respond  ing  to  the  particular 

wave-mode  coordinate,  e^_ .  is  a  unit  vector  pointing  in  the  direction  of 

t 

increasing  x.,  where  x^  is  the  (local)  x-axls  of  the  i  member  attached 
to  the  joint,  and  is  the  outward  unit  normal  to  the  joint  at  the  point 

where  the  member  is  attached  to  the  joint.  Eqn.  (20)  will  be 

satisfied  if  the  x^-axis  points  "into"  the  joint  and  the  wave-mode 
coordinate  propagates  in  the  direction  of  increasing  x^,  or  if  the 
x  .-axis  points  "away  from"  the  joint  and  the  wave-mode  coordinate  propa¬ 
gates  in  the  direction  of  decreasing  x^.  A  wave-mode  coordinate  in  the 
itn  member  attached  to  the  joint  represents  a  wave  leaving  the  joint  if 


points  "into"  the  joint 


Eqn.  (21)  will  be  satisfied  if  the  x^-axis 
and  the  wave-mode  coordinate  propagates  in  the  direction  of  decreasing 
x^,  or  if  the  x^-axis  points  "away  from"  the  joint  and  the  wave-mode 
coordinate  propagates  in  the  direction  of  increasing  x^.  Fig.  1 
illustrates  schematically  wave-mode  coordinates  entering  and  leaving 
a  generic  joint. 

The  scalar  equations  in  eqn.  (18)  may  now  be  rearranged  to  give 

w 

—out 
D 

w  . 

—  m 

wTiere  w  is  a  vector  containing  all  the  wave-mode  coordinates  leaving 
—out 

the  joint,  w  is  a  vector  containing  all  the  wave-mode  coordinates 
entering  the  joint,  and  D  is  the  matrix  obtained  by  rearranging  the 
scalar  equations  of  eqn.  (18)  into  the  form  of  eqn.  (22).  (The  matrix 
manipulations  leading  to  eqn.  (22),  and  all  the  other  matrix  manipu¬ 
lations  described  in  this  section,  will  become  clearer  in  the  examples 
given  below.)  Eqn.  (22)  can  be  written  as 


D  w  +  D .  w .  =  f 

—out  —out  —  in  —in  —ext 


(23) 


where  D  is  the  submatrix  of  D  which  multiplies  w  .  and  D.  is  the 
—out  —  —out  —in 

submatrix  of  JD  which  multiplies  .  Premultiplication  of  both  sides 
of  eqn.  (23)  by  (assuming  that  exists)  gives 

w  +  D  ^  D .  w.  =  D  ^  f 
—out  —out— in— in  —out  —  ext 


(24) 


The  matrix  is  called  the  scattering  matrix  of  the  joint,  and  the 

matrix  G(ijj)  is  called  the  wave-mode  generation  matrix  of  the  joint. 

Eqn .  (25),  which  is  merely  a  manipulated  form  of  eqn.  (17),  is  the 
final  form  of  the  description  of  the  dynamics  of  a  joint  in  terms  of 
wave-mode  coordinates.  Eqn.  (25)  may  be  considered  as  an  input-output 
relationship  for  the  joint,  and  states  that  outgoing  waves  may  be 
generated  by  incoming  waves  via  the  scattering  matrix  j>(oo) ,  or  may  be 
generated  by  external  forces  (or  moments)  via  the  wave-mode  generation 
matrix  jQ(to)  .  The  scattering  matrix  S(u>)  contains  the  transmission  and 
reflection  coefficients  which  the  wave-mode  coordinates  encounter  as 
they  enter  the  joint.  If  the  transmission  and  reflection  coefficients 
in  S^(uj)  depend  on  the  frequency  co,  the  scattering  at  the  joint  is  called 
dispersive.  If  the  transmission  and  reflection  coefficients  in  S>(w) 
are  independent  of  frequency,  the  scattering  is  called  nondispersive . 

With  a  scattering  matrix  and  a  wave-mode  generation  matrix  for 
each  joint  in  a  lattice,  and  a  wave-mode  matrix  and  a  wave-mode  propagation 
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matrix  for  each  member  in  a  lattice,  it  is  possible  in  principle  to 
analyze  wave  propagation  in  lattice  structures  of  arbitrary  complexity. 
Following  several  intermediate  examples,  a  complete  simple  example  of 
such  an  analysis  is  given  in  the  next  section. 


EXAMPLES 
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where  the  transfer  matrix  T^  is  given  by  [3] 


cos6 


T 


sin6 

0 


j  sin  6 
EA  0 


COS0 


where 


Z  = 


-  x 


1 


y  =  pa 


The  eigenvalues  of  the  matrix  T  in  eqn.  (31)  are  [8] 


=  cos0 


isin0 
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e 


=  cos0  +  isin0 


ie 

e 


Therefore,  the  propagation  constants  given  by  eqn.  (12)  are 


where 


(45) 


The  wave-mode  vector  w^  at  x  =  X2  and  the  wave-mode  vector  w^  at 
x  =  x^  are  related  according  to  eqn.  (8)  by 


W2  =  W(u')w^ 


(46) 


Writing  out  eqn.  (46)  gives 


0 1 

r-  +  1 

~  "f  ‘ 

w„ 

W., 

2 

= 

l 

w2 

0 

w 

1 

(47) 


Eqn.  (47)  is  the  diagonalized  form  of  the  transfer  matrix  relationship 

given  by  eqn.  (30) .  The  components  w^  and  w*  of  the  wave-mode  vectors 

W2  and  w^  in  eqn.  (47)  are  written  with  a  plus  (+)  superscript  since  the 

—  -f  ■—  -f 

propagation  constant  -y  which  corresponds  to  and  w^  has  a  negative 
imaginary  part,  indicating  that  and  w^"  represent  waves  which 

travel  in  the  direction  of  increasing  x.  Similarly,  the  components 

and  w^  are  written  with  a  minus  (-)  superscript  since  the  propa¬ 
gation  constant  y  which  corresponds  to  w  and  w 


has  a  positive 


imaginary  part,  indicating  that  w^  and  represent  waves  which 
travel  in  the  direction  of  decreasing  x. 

From  eqns .  (4)  and  (43),  the  wave-mode  coordinates  for  the  longi¬ 
tudinal  rod  are  given  in  terms  of  the  physical  components  of  the  state 
vector  by 


u(x,aO 

F(x,io) 


(4b) 
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u(x,tj)  =  W+(x,(jL>)  +  w  (x,u) 

F(x,co)  =  -iRu.'  w +  (x,c)  +  iRu;  w  (x,iu) 


(52) 

(53) 


Eqn.  (52)  expresses  the  familiar  fact  that  the  displacement  in  an 
elastic  longitudinal  rod  may  be  written  as  the  sum  of  a  wave  travelling 
in  the  direction  of  increasing  x  and  a  wave  travelling  in  the  direction 
of  decreasing  x.  Eqn.  (53)  states  that  the  force  (or  stress)  in  the  rod 
can  also  be  written  as  the  sum  of  the  same  two  waves,  scaled  (in  the 
frequency  domain)  by  the  appropriate  factors. 

The  eigenvectors  given  by  eqns .  (39)  and  (40)  remain  eigenvectors 
when  multiplied  by  any  constant.  By  choosing  the  first  components  of 
the  eigenvectors  to  be  unity,  the  wave-mode  coordinates  in  eqn.  (48) 
are  given  the  dimensions  of  displacement.  If  the  second  components 
of  the  eigenvectors  are  chosen  to  be  unity,  the  resulting  wave-mode 
coordinates  will  have  the  dimensions  of  force.  Both  sets  of  wave-mode 
coordinates  will  satisfy  eqn.  (47). 

The  propagation  constants  given  by  eqns.  (37)  and  (38)  are  purely 
imaginary.  Thus,  waves  on  an  elastic  longitudinal  rod  are  not 
attenuated  as  they  propagate.  Further,  the  imaginary  parts  of  the 
propagation  constants  given  by  eqns.  (37)  and  (38)  are  linearly  pro¬ 
portional  to  frequency,  and  thus  waves  in  an  elastic  rod  are  nondis- 
persive.  The  phase  velocity  of  the  waves,  which  as  stated  above  is 
equal  to  the  inverse  of  the  constant  of  proportionality,  is,  as  is 
well  known,  /E/p  . 
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L„  imple  2 :  Joint  Coupling  Matrix  for  a  Rigid  One-Dimensional 
Joint  with  Mass 

In  this  example,  the  joint  coupling  matrix  for  a  one-dimensional 
rigid  joint  is  derived.  A  rigid  one-dimensional  joint  with  mass  m  is 
shown  in  Fig.  3a  Two  members,  denoted  by  1  and  2,  are  attached  to 
the  joint.  The  joint  is  called  one-dimensional  because  all  members 
attached  to  the  joint  lie  along  a  line,  and  the  joint  is  constrained 
to  move  only  along  that  line.  The  joint  is  subjected  to  an  external 

force  F  .  Fig.  3a  also  shows  local  coordinate  axes  x^  and  for 

members  1  and  2,  respectively,  and  a  joint  coordinate  axis  x^.  Fig.  3b 

shows  the  components  u^  and  (i  =  2,3)  of  the  state  vectors  of  the  two 

members  at  their  respective  points  of  contact  with  the  joint,  and  the 
components  and  F^  of  the  state  vector  of  the  joint. 

Since  the  joint  is  rigid,  the  geometric  compatibility  requirement 
at  the  joint  is 

u2  =  u3  =  Uj  (54) 

Therefore, 

G2  =  u3  =  ua  (55) 

from  which  it  follows  that 

-u2  +  u3  =  0  (56) 

Applying  Newton's  Second  Law  to  the  joint  gives,  in  terms  of  the 
Fourier  transforms  of  the  displacements  and  forces, 

-F2  +  F^  +  Fj  =  -raui^Uj  (57) 
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Substitution  of  eqn.  (54)  into  eqn.  (57)  gives 


-F^  +  ¥  j  +  F  j  =  -mu)2u2 


(58) 


or 


-mw2a2  +  F2  -  F3  =  FJ 


(59) 


Eqns.  (56)  and  (59)  can  be  written  as 

r 


G2 

-1  0  1  O'"1 

0 

2  I  < 

_ 

►  — 

-mu  1  0  -1 

u„ 

F 

1—  —1 

3 

J 

F. 

3 

k  * 

(60) 


Eqn.  (60)  is  in  the  form  of  the  general  joint  coupling  relationship 
given  by  eqn.  (15),  with 


B  = 


-1 


0  1  O' 


-mu  1  0  -1 


(61) 


Eqn.  (61)  gives  the  joint  coupling  matrix  for  the  one-dimensional 
joint  of  Fig.  3. 

This  simple  example  illustrates  the  general  procedure  for  deriving 
a  joint  coupling  matrix.  The  joint  coupling  relationship  for  rigid 
joints  consists  essentially  of  geometric  compatibility  requirements 
and  dynamic  requirements,  written  in  terms  of  Fourier  transforms  and 
in  matrix  form.  In  [7],  the  procedure  used  to  derive  the  joint  coupling 
matrix  of  eqn.  (61)  is  applied  to  more  complex  two  and  three-dimensional 
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rigid  joints. 

Example  3:  Scattering  Matrix  for  the  Junction  of  Two  Elastic 
Longitudinal  Rods 

Two  elastic  longitudinal  rods  connected  by  a  rigid  one-dimensional 
joint  with  mass  m  are  shown  in  Fig.  4a.  Rod  1  has  mass  density  p  , 
cross-sectional  area  A^,  and  elastic  modulus  E^.  Rod  2  has  mass  density 
p 2>  cross-sectional  area  ,  and  elastic  modulus  ■ 

As  discussed  in  the  previous  example,  the  dynamics  of  the  joint  in 
Fig.  4a  are  described  by  the  joint  coupling  relationship  given  by 

eqn.  (60).  In  this  example,  the  process  of  transforming  eqn.  (60)  into 

a  scattering  relationship  of  the  form  of  eqn.  (25)  is  illustrated. 

From  eqn.  (51),  the  state  vectors  at  points  2  and  3  in  Fig.  4 

are  related  to  the  wave-mode  coordinates  at  points  2  and  3  by 


U2 

F_ 

2 

-iR^oj  iR^w 


'  — 

U3 

?3 

1  1  l»3+' 


—  i  R  1R2w  [w 


where 


R1  ~  A1  v  plEl 


^2  A2  v 
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The  waves  represented  by  w^"  ,  ,  w^+  ,  and  w^  are  shown  schematically 


by  arrows  in  Fig.  4b.  Substitution  of  eqns.  (62)  and  (63)  into  eqn . 
(60)  gives 


0  1 


-mu  1  0 


0 

-1 
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r  1 
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-  +1 
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iR^uo 
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W3 
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iR2w 

■ 
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(66J 


which  can  written  as 


-1  0  10 
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-rrw  10-1 
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w2 

-iR^u' 

iR^iu 

0 

0 

*2_ 

-  + 
w 

3 

0 

0 

0 
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►  = 
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-R2u) 

iR2U) 

w 

lw3  J 

of  eqn. 

(18). 

Multiplying 

out  eqn. 

(67) 

(67) 


rw  + 

2 


-1 


,  2  2 

(-mu  -iR^u)  (-mu  +iR^u)  iR^t-)  -iR^' 


w„ 


>  = 


-  + 
W3 


(68) 


^  w^  J 


-  + 


Note  from  Fig.  4b  that  w^  and  w^  represent  waves  leaving  the  joint. 


-  + 


and  w0  and  w^  represent  waves  entering  the  joint.  Rearranging  the 
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m  w  ^  l%  .v.  ->r.  vr  v  < 


yr L»  J.  W  M  w 


!• 
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order  of  the  scalar  equations  contained  in  eqn.  (68)  gives 

--1  1 

2 

(-mm  +iR^w)  iR^a 

Eqn.  (69)  is  in  the  form  of  eqn.  (22)  with 


1  1 


(-me  -iR  a) 


-iR„ 


-  + 
W3 
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w2 


v.w3  s 
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and 


w 

—out 


Eqn.  (69)  is  eqivalent  to 


-1  1 

*2 

2 

-  + 

(-mw  +iR1w)  iR?oj 

Eqn.  (72)  is  in  the  form  of  eqn.  (23),  with 


D  _ 
— out 


-1 


(-mu  +  iR^oj)  iR2u 


and 


D . 
—in 


-1 


(-mu  -  iRjUi)  -iR^w 


The  inverse  of  the  matrix  I)  given  by  eqn.  (73)  is 

out 


-iR2^ 


D  1 
—out 
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-mu  -t-iR^oj+iR^ 


2 

-mu)  +iR^W 


-mu  +  iR^uH-iR^u 


-mu  +iR^(y+iR2w 


-mu  -KR^on-iR^ 


Premultiplication  of  both  sides  of  eqn.  (72)  by  DQUt;  gives,  a 
some  manipulation, 
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Eqn .  (76)  is  in  the  form  of  eqn.  (25),  with  the  scattering  matrix  j3((0 


given  by 


SCO  = 


mix.  -t-iR^U'-iR^ 


2iR2u) 


-me.'  +iR^0J+iR2u 


-mu>  +iR^u)+iR2to 


2iR  u 


mu,  -iR  uH-iR2u) 


-mu'  +iR^ui+iR2u 


-mu  +iR^uH-iR2ui 


and  the  wave-mode  generation  matrix  C(u)  given  by 


-mu  +iR^u*+iR2o) 


-mu  “+i  R^u+iR2ui 


G(u>)  = 


-mu.  +iR^w 


-mu  +iR^uH-iR2Uj 


-mu  +iR^u'+iR2w 
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Eqn.  (76)  is  the  description  of  the  dynamics  of  the  joint  of 
Fig.  4  in  terms  of  wave-mode  coordinates.  Note  that  the  description  of 
the  joint  dynamics  in  terms  of  the  physical  components  of  the  state 
vectors  (eqn.  (60))  contains  only  information  about  the  dynamics  of  the 
joint.  On  the  other  hand,  the  description  of  the  joint  dynamics  in 
terms  of  wave  coordinates  given  by  eqn.  (76),  which  implicitly  contains 
the  joint  coupling  matrix  given  by  eqn.  (61)  and  the  wave-mode  matrices 
given  by  eqns.  (62)  and  (63),  contains  information  about  both  the  joint 
dynamics  and  the  dynamics  of  the  members  attached  to  the  joint. 

The  quantity 

2  . 

mui  +iR^oj-iR2U} 

2 

-ma.'  +iR^'jJ+iR2<jJ 

in  the  scattering  matrix  S^(oj)  given  by  eqn.  (77)  is  the  reflection 
coefficient  for  w2+ »  which  represents  a  wave  incident  upon  the  joint 
from  the  left,  and  the  quantity 

2iR^io 

2 

-mto  +iR^oj+iR2W 

is  the  transmission  coefficient  for  W2+  .  Similarly,  the  quantity 
2 

mcu  -iR  w+i^u) 

2 

-mai  +iR^uj+iR2W 

is  the  reflection  coefficient  for  w^  ,  which  represents  a  wave  incident 
upon  the  joint  from  the  right,  and  the  quantity 
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-mw  +iR^aH-iR20J 


is  the  transmission  coefficient  for  .  Note  that  since  each  wave-mode 
coordinate  contains  information  about  both  displacement  and  force  (or 
stress) ,  it  is  not  necessary  when  using  wave-mode  coordinates  to  derive 
one  set  of  reflection  and  transmission  coefficients  for  displacement 
waves  and  a  different  set  of  reflection  and  transmission  coefficients 
for  force  (or  stress)  waves,  as  is  done,  for  example,  in  [9].  Since  the 
wave-mode  coordinates  used  here,  which  are  defined  by  eqn.  (48),  have 
the  dimensions  of  displacement,  the  reflection  and  transmission 
coefficients  in  the  scattering  matrix  j[(to)  given  by  eqn.  (77)  correspond 
to  the  displacement  reflection  and  transmission  coefficients  given  in  [9] 
The  transmission  and  reflection  coefficients  contained  in  the 
scattering  matrix  _S (oo)  given  by  eqn.  (77)  depend  on  the  frequency  0). 
Therefore,  this  is  an  example  of  dispersive  scattering.  As  the  frequency 
a)  becomes  very  large,  the  scattering  matrix  j5(u))  given  by  eqn.  (77) 
approaches 


S  (to) 


(79) 


Thus,  for  very  high  frequencies,  the  scattering  is  nondispersive ,  and 
the  reflection  coefficients  are  equal  to  -1,  while  the  transmission 
coefficients  are  equal  to  zero.  Thus  for  very  high  frequencies,  the 


joint  in  Fig.  4  behaves  like  a  rigid  boundary. 


If  the  mass  of  the  joint  in  Fig.  4  is  zero,  the  scattering  matrix 
j>(u)  given  by  eqn.  (77)  becomes 


S(w)  = 


R1  ~  R2 
Rl  +  R2 


™1 

R,  +  R, 


2R2  ~[ 
+  R- 


^2 _ ^1 

Rl  +  R2 


The  scattering  matrix  given  by  eqn.  (80)  contains  the  commonly  derived 
(see,  for  example,  [9])  reflection  and  transmission  coefficients  for 
displacement  waves  at  the  junction  of  two  collinear  longitudinal  rods. 


Example  4:  Analysis  of  a  One-Dimensional  Lattice  Structure 


In  this  fourth  and  final  example,  which  utilizes  the  results  of 
the  previous  three  examples,  the  one-dimensional  lattice  of  Fig.  5  is 
analyzed  using  wave-mode  coordinates  and  scattering  matrices.  The 
lattice  consists  of  two  elastic  longitudinal  rods  connected  by  a  rigid 
one -dimensional  joint.  It  is  assumed  that  the  joint  is  massless.  Rod 
1  has  mass  density  p^,  cross-sectional  area  A^,  elastic  modulus  E^,  and 
length  Rod  2  has  mass  density  p^>  cross-sectional  area  A9,  elastic 

modulus  E^,  and  length  The  joint  is  subjected  to  a  given  force 

,  and  it  is  desired  to  find  the  resulting  force  F^(t)  at  point  1. 
The  scattering  matrix  relationship  at  the  joint  is  given  by  eqn. 
(76)  as 
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ia)(R1+R2) 


where 


Ri 
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+  R2 

2R1 

(82) 


(83) 


(84) 


R1  +  R2 


(85) 


and  R  and  R0  are  given  by  eqns.  (64)  and  (65),  respectively. 


The  components  of  the  wave-mode  vector  w^  at  point  1  in  Fig.  5 


satisfy  the  relationship 


(w^  )  =  ro(w~  ) 


(86) 


where  r  is  the  reflection  coefficient  at  the  left-hand  boundary  (point 


1)  of  rod  1.  The  quantity  r  may  be  derived  from  the  boundary  conditions 


at  point  1.  If,  for  example,  point  1  is  a  fixed  end,  r  =  -1,  and  if 


point  1  is  a  free  end,  r  =1  [10],  Similarly,  the  components  of  the 

o 


wave-mode  vector  w^  at  point  4  in  Fig.  5  satisfy  the  relationship 


(w4  )  =  r3^4+  > 


(87) 
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where  is  the  reflection  coefficient  at  the  right-hand  boundary 
(point  4)  of  rod  2.  Again,  the  quantity  r^  may  be  derived  from  the 
boundary  condition  at  point  4;  if  point  4  is  a  fixed  end,  r^  =  -1, 


and 


if  point  4  is  a  free  end,  r^  =  1.  The  quantities  r^  in  eqn.  (86)  and 
r^  in  eqn.  (87)  may  be  considered  as  1  x  1  scattering  matrices.  It  is 
assumed  that  the  quantities  r^  and  r^  are  independent  of  frequency. 

Eqns .  (81),  (86)  and  (87)  are  the  scattering  relationships  for  the 
lattice  of  Fig.  5. 

The  components  of  the  wave-mode  vector  w^  at  point  2  in  Fig.  5  and 
the  components  of  the  wave-mode  vector  w^  at  point  1  in  Fig.  5  are 
related  by  a  wave-mode  propagation  matrix  according  to  eqn.  (47)  ar 
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where 


Y1  =  1U\/^ 


(89) 
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where 


Y2  =  Xl\Z 


Eqns.  (88)  and  (90)  are  the  wave-mode  propagation  relationships  for  the 
lattice  of  Fig.  5. 

Eqns.  (81),  (86),  (87),  (88)  and  (90)  contain  eight  scalar 

equations  for  the  eight  unknowns  w*  ,  w^  ,  w^+  ,  w^  ,  ,  w^  s  w^ 

and  w.  .  These  equations  mav  be  solved  for  w+  and  w  to  give  [ 11 ] 
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Eqns.  (9~)  and  (93)  give  the  wave-mode  coordinates  w^  and  w^  in  terms 
of#.  The  physical  components  of  the  state  vector  at  point  1  are 
given  in  terms  of  the  wave-mode  coordinates  at  point  1  by  eqn.  (51)  as 


-iR,u'  iR,e>  (  w, 
i  1  _j  '  1 


Substitution  of  eqns.  (92)  and  (93)  into  the  second  scalar  equation  of 
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eqn.  (94)  gives 
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Eqn.  (95)  gives  the  Fourier  transform  of  F^(t)  in  terms  of  the  Fourier 
transform  of  t)  •  Thus,  the  quantity  in  the  square  brackets  in 

eqn.  (95)  is  the  transfer  function  between  the  force  ^  and  the  force 
F^.  As  will  be  shown  below,  the  transfer  function  in  eqn.  (95)  can  be 
written  as  an  infinite  series  of  constants  multiplied  by  pure  delays. 
It  is  noted  in  passing  that  the  values  of  to  for  which  the  denominator 
of  the  transfer  function  in  eqn.  (95)  vanishes  are  the  natural 
frequencies  of  the  lattice  of  Fig.  5,  and  that  there  is  an  infinite 
number  of  such  frequencies.  Thus,  the  lattice  of  Fig.  5  is,  in  the 
language  of  linear  system  theory,  an  "infinite-order  system." 

In  order  to  find  F^(t),  it  is  necessary  first  to  know  ,  and 

then  to  perform  an  inverse  Fourier  transform  on  eqn.  (95) .  In  this 
example  it  is  assumed  that 


where  6(t)  is  the  Dirac  delta  function,  so  that  [2] 
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eqn.  (99)  can  be  written  as  [11] 
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Now,  recalling  the  definitions  of  y  and  Y2  given  by  eqns.  (89)  and  (91), 

i-  _  i^T  .  .  .  „ , 

and  the  fact  that  the  inverse  Fourier  transform  of  e  is  o(t-T)  |2J, 

the  inverse  Fourier  transform  of  eqn.  (101)  can  be  taken  term  by  term 


to  give 
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where 
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and  A(l)  is  a  time-shift  factor  defined  by 


f ( t) A (t)  =  f(t-T) 
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Writing  out  the  first  few  terms  of  eqn.  (130)  gives 
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Eqn.  (103),  which  is  an  infinite  series  of  impulses  scaled  in  amplitude 
by  reflection  and  transmission  coefficients  and  delayed  by  multiples 
of  the  times  T ^  and  is  the  desired  solution  for  F^(t). 

Each  of  the  impulses  in  ean.  (107)  can  be  interpreted  physically 

by  considering  the  propagation  of  the  initial  impulse  applied  to  the 

joint.  The  initial  impulse  &  causes  an  impulse  to  propagate  to  the 

left  towards  point  1,  and  also  causes  an  impulse  to  propagate  to  the 

right  towards  point  A.  After  a  time  delay  i  ^ ,  the  impulse  which 

initially  propagates  to  the  left  arrives  at  point  1,  where  it  is 

reflected  and  scaled  in  amplitude  bv  the  reflection  coefficient  r  . 

o 

The  reflected  impulse  arrives  at  the  joint  after  an  additional  time 
delay  T  ^ ,  and  at  the  joint  it  is  partially  reflected  back  towards  point 
1,  and  partially  transmitted  towards  point  A.  The  reflected  impulse  is 
scaled  in  amplitude  by  the  reflection  coefficient  r^,  and  the  trans¬ 
mitted  impulse  is  scaled  in  amplitude  by  the  transmission  coefficient 
t7.  The  impulse  which  initially  propagates  from  the  joint  to  the  right 
arrives,  after  a  time  delay  i „ ,  at  point  A,  where  it  is  reflected  and 
scaled  in  amplitude  by  the  reflection  coefficient  r^.  The  reflected 
impulse  arrives  at  the  joint  after  an  additional  time  delay  X7,  and  at 
the  joint  it  is  partially  reflected  back  tow-ards  point  A,  and  partially 
transmitted  towards  point  1.  The  reflected  impulse  is  scaled  in  ampli¬ 
tude  by  the  reflection  coefficient  r,,  and  the  transmitted  impulse  is 
scaled  in  amplitude  by  the  transmission  coefficient  t^.  This  process 
of  propagation,  reflection,  and  transmission  continues  indef  initelv . 
Since  the  propagation  constants  and  the  scattering  matrices  for  the 
lattice  of  Fig.  5  are  nond ispersive,  the  initial  impulse  applied  to  the 
joint  maintains  its  identitv  as  an  impulse  as  it  propagates  and  as  it 


is  reflected  and  transmitted. 

The  first  infinite  series  in  eqn.  (107)  represents  those  impulses 
which  initially  propagate  from  the  joint  to  the  left,  and  then  remain 
in  rod  1  due  to  successive  reflections  at  point  1  and  at  the  joint. 

The  first  product  of  two  infinite  series  in  eqn.  (107)  represents 
those  impulses  which  initially  propagate  from  the  joint  to  the  right, 
and  then,  after  an  arbitrary  number  of  reflections  in  rod  2,  are  trans¬ 
mitted  through  the  joint  into  rod  1,  and  remain  thereafter  in  rod  1. 

The  second  product  of  two  infinite  series  in  eqn.  (107)  represents 
those  impulses  which  initially  propagate  from  the  joint  to  the  left, 
reflect  an  arbitrary  number  of  times  in  rod  1,  are  transmitted  through 
the  joint  into  rod  2,  reflect  an  arbitrary  number  of  times  in  rod  2, 
are  transmitted  through  the  joint  into  rod  1,  and  remain  thereafter 
in  rod  1.  The  third  product  of  two  infinite  series  in  eqn.  (107) 
represent  those  impulses  which  initially  propagate  from  the  joint  to 
the  right,  reflect  an  arbitrary  number  of  times  in  rod  2,  are  trans¬ 
mitted  through  the  joint  into  rod  1,  reflect  an  arbitrary  number  of 
times  in  rod  1,  are  transmitted  through  the  joint  into  rod  2,  reflect 
an  arbitrary  number  of  times  in  rod  2,  are  transmitted  through  the  joint 
into  rod  1,  and  remain  thereafter  in  rod  1.  The  remaining  terms  in 
eqn.  (103)  can  be  interpreted  in  a  similar  manner. 

The  response  F^(t)  given  by  eqn.  (103)  is  obtained  here  by  using 
wave-mode  coordinates  and  scattering  matrices  in  the  frequency  domain 
to  obtain  an  expression  F^(w),  and  then  performing  an  inverse  Fourier 
transform  to  obtain  F^(t) .  The  response  F^(t)  can  also  be  obtained  by 
a  direct  consideration  of  propagation,  reflection,  and  transmission  in 


the  time  domain.  Such  an  approach  is  taken  in  [13].  The  response 
F^(t)  given  here  agrees  exactly  with  the  response  which  is  obtained 
using  the  methods  given  in  [13] .  In  fact,  the  numerical  coefficients 
which  appear  in  the  various  infinite  series  of  eqn.  (103),  and  which 
occur  here  naturally  as  a  part  of  the  process  of  inverse  Fourier 
transformation,  account  for  the  existence  of  "equivalent  paths"  from 
the  input  location  (the  joint)  to  the  response  location  (point  1).  The 
concept  of  equivalent  paths  in  a  lattice  structure  is  discussed 
thoroughly  in  Appendix  A  of  [13],  and  some  sketches  of  equivalent 
paths  in  the  lattice  of  Fig.  5  are  given  in  [11].  The  advantage  of 
the  approach  taken  here,  which  uses  wave-mode  coordinates  and  scattering 
matrices,  is  that  this  approach  may  be  extended,  using  the  same  formal 
procedures,  to  the  analysis  of  two  and  three-dimensional  lattice 


structures. 


CONCLUSIONS 


The  examples  given  in  the  previous  section  illustrate,  in  principle 
all  the  steps  required  for  an  analysis  of  wave  propagation  in  an 
arbitrary  lattice  structure.  The  basic  ingredients  of  such  an  analysis 
are  a  wave-mode  matrix  Y(u,)  and  a  wave-mode  propagation  matrix  W(to) 
for  each  member  of  the  lattice  and  a  scattering  matrix  S^(oj)  and  a 
wave-mode  generation  matrix  C(u})  for  each  joint  in  the  lattice.  A 
wave-mode  propagation  relationship  such  as  eqn.  (7)  applied  to  each 
member  and  a  scattering  relationship  such  as  eqn.  (25)  applied  to  each 
joint  provide  a  set  of  equations  which  can  be  solved  for  the  wave-mode 
vector  w  at  any  point  in  the  lattice.  The  state  vector  at  any  point 
in  the  lattice  can  then  be  obtained  through  the  use  of  the  wave-mode 
transformation  given  by  eqn.  (5).  Finally,  individual  physical  com¬ 
ponents  of  the  state  vector  jj  may  be  obtained  through  the  process  of 
inverse  Fourier  transformation. 

The  major  difficulty  in  the  procedure  just  described  is  the 
evaluation  of  the  inverse  Fourier  transform.  Efficient  numerical 
techniques  for  evaluating  complicated  inverse  Fourier  transforms  are 
necessary  if  the  methods  described  in  this  paper  are  to  be  used 
successfully  in  the  analysis  of  large  and  complex  lattice  structures. 
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Two  elastic  longitudinal  rods  attached  to  a  rigid 
one-d imensional  joint. 


